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MAGNETOROTATIONALLY DRIVEN SUPERNOVAE AS THE ORIGIN OF 
EARLY GALAXY i?-PROGESS ELEMENTS? 
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ABSTRACT 

We examine magnetorotationally driven supernovae as sources of r-process elements in the early 
Galaxy. On the basis of thermodynamic histories of tracer particles from a three-dimensional magne- 
tohydrodynamical core-collapse supernova model with approximated neutrino transport, we perform 
nucleosynthesis calculations with and without considering the effects of neutrino absorption reactions 
on the electron fraction [Y^) during post-processing. We find that the peak distribution of Y^. in the 
ejecta is shifted from ~ 0.15 to '^ 0.17 and broadened toward higher Yg due to neutrino absorption. 
Nevertheless, in both cases the second and third peaks of the solar r-process element distribution can 
be well reproduced. The rare progenitor configuration that was used here, characterized by a high 
rotation rate and a large magnetic field necessary for the formation of bipolar jets, could naturally 
provide a site for the strong r-process in agreement with observations of the early galactic chemical 
evolution. 

Subject headings: magnetohydrodynamics (MHD) — nuclear reactions, nucleosynthesis, abundances 
— stars: magnetic field — stars: neutron — stars: rotation — supernovae: general 



1. INTRODUCTION 

The supernova (SN) mechanism explaining how mas- 
sive stars end in a central collapse to a neutron star (or 
a black hole) and explosive ejection of the outer layers 
is still debated. It has been related to neutrino e mission 
from the hot collapsed core and accreted matter ^Beth3 
[1990), and this (regular) scenario is again becoming 
more and more promising with recent multi-D hydro- 
dynamic approac hes and im proved ( spectral) neutrino 
transport (e.g.. iMiiller et a l. 201 0; iLentz et al.l 120121: 
iLiebendorfer et all l201cinBlirrows et all 120111) . The 
equation of state (EoS) caused explosion scenario has 
found some recent revival, based on the quark-hadron 
phase transition at supernuclear densities (jSagert et al.l 
l2009f l. but its working depends on specific choices of 
EoS properties in a very narrow parameter range. Other 
options are related to rotation and magnetic fields, a 
topic which has been discussed for more than 30 years, 
but required three-dimensional (3D) modeling and had 
therefore only been addressed with limited success in the 
early days. The major outcome was that high rotation 
rates and (possibly unrealistically high) magnetic fields 
were required to launch explosions. The question is 
whether such magnetic fields can be attained during 
collapse with rotation and on wh i ch timescale after 
collapse (lleBlanc fc Wilsonl 119701: iMeier et all 119761: 



Bisnovatvi-Kogan et al.l 1 19761 IMiiller fc Hillebrandd 
19791 and references th erein). This topic has recently 
been re-addressed by iMoiseenko fc Bisnovatvi-KoganI 
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(|2008[ ) and iTakiwaki et"al] ()2009[ ). In the present 
Letter we want to present results obta ined with our 
3D m agneto-hydrodynamics code FISH (jKappeli et al.l 

unni). 

There is another reason why such magnetohydrody- 
namically driven explosions are of possible interest: the 
search for the site of a (full) r-process early in galac- 
tic evolution. The initially postulated neutrino wind 
in re g ular (neutrino-dr i ven) supernovae dHoffman et al.l 
[19971 : IMever fc Bro^ [19971 : iFarouai et all [20l0h Ca^ 
did not result in hav i ng suffic i ently high entropies 
([Thompson et all 120011: iWanaiol 120061: lArcones et all 
120071 ). and (b) in addition, the innermost e jecta turned 
out t o be proton-rich instea d of neutron- rich ([Pruet et al.l 
120061 : iFrohlich et al.l l2006f ) , a situation w hich has now 
also been shown in lon g-term simulations ([Fischer et al.l 
120101 : IHiidepohl et al.l 12010). Electron-capture super- 
novae, which explode without a long phase of accretion 
onto the proto-neutron sta r (PNS), apparently provide 
more favorable conditions ([Wanaio et al.l [201 ID . How- 
ever, also the proton to nucleon ratios obtained under 
such conditions do not support a strong r-process, which 
successfully reproduces the platinum peak of r-elements 
around A — 195. A similar result was recently obtained 
for core-collapse supernova (CCSN) explosions triggered 
by a quark-hadron phase transition during the early post- 
bounce ph ase when investigating their detailed nucle- 
osynthesis ([Nishimura et al.ll2011l ). Both types of events 
might support a weak but not a full r-process. 

Neutron star mergers have been shown to be power- 
ful sources of r-process matter, in fact ejecting a factor 
of 100 to 1000 more r-process material than required on 
average from CCSNe, if those would have to explain so- 
lar r-process element abundances. This would actually 
support the large scatter of Eu/Fe found in very metal- 
poor stars. The only problem is that it might be hard 
to explain the early appearance of r-proce ss matter for 
metallicities at and below [Fe/H] — 3 ([Argast et al.l 



Winteler et al. 



120041) . Some recent studies, which include the fact that 
our Galaxy is possibly the result from smaller merg- 
ing subsystems (with different star formation rates) have 
been expected to show a way out of this dilemma. If 
this cannot be solved, we need another strong r-process 
source already at low metallicities, and possibly jets from 
rotating core coll apses with strong magnetic fields could 
be the solution |Cameronl 120031 : iNishimura et al.l [20061 : 
iFuiimoto et al.ll2008[) . 

The present Letter has the aim to explore the results 
from our 3D magnetohydrodynamic (MHD) calculations, 
which lead to bipolar jet ejection. The following Sec- 
tion [5] will discuss the initial models and the explosion 
dynamics; Section |3] will present nucleosynthesis results. 
Section |4] is devoted to a discussion of uncertainties and 
an outlook on future investigations. 

2. 3D MHD-CCSN MODEL 

The calculation presented here was performed with the 
comp utational setup sirnilar to our previous investiga - 
tions (|Liebend6rfer et al.l 120051 : iScheidegger et al.ll2010f ). 
The initially innermost (600 km)'^ of the massive star are 
covered by a 3D Cartesian domain uniformly discretized 
by 600^ cells, resulting in a 1 km resolution, that is em- 
bedded in a spherically symmetric domain encompassing 
the iron core and parts of the silicon sh ell. The mag- 
netic fluid is evolved with the FISH code (|Kappeli et al.l 
1201 It ), solving the ideal MHD equations. The spheri- 
cally symmetric domain is evolved with the AGILE code 
(|l7iebendorfcr et al. 2002) . The gravitational potential is 
approximated by an effective axisymmetric mass distri- 
butio n that includes gen eral relativist ic monopole correc- 
tions (iMarek et al.ll2006f ) . We use the lLattimer fc Swestvl 
(|1991h EoS with nuclear compressibility 180 MeV. We 
have included a Lagrangian component in the form of 
tracer particles which are passively advected with the 
flow. They record the thermodynamic conditions of a 
particular fluid element and serve as input to the post- 
processing nucleosynthesis calculations. 

The transport of the electron neutrinos and anti- 
neutrinos is approximated by a 3D spectral leak- 
age scheme, based on previou s gray leakage schemes 
(jRosswog fc Liebendorfeii [20031 and references therein). 
The neutrino energy is discretized with 12 geometri- 
cally increasing energy groups spanning the range E,^ = 
3 200 MeV. The amount of energy and particles lo- 
cally released is calculated for each bin as an interpola- 
tion between the diffusive rates and the (free streaming) 
production rates, depending on the local neutrino optical 
depth. For the computation of the spectral optical depth 
we have used a ray-by-ray axisymmetric approximation, 
calculated on a polar grid encompassing the full 3D 
Cartesian domain discretized uniformly with 1km radial 
spacing and 30 angular rays covering the full [0, n] realm. 
All fundamental neutrino reactions have been included 
(neutrino scattering on nucleons and nuclei, neutrino ab- 
sorption/emission on nucleons and nuclei) , providing de - 
tailed spectral emissivities and opacities (|Bruennlll985l ). 
Inside the neutrinosphere, weak equilibrium is assumed 
and trapped neutrinos are modeled accordingly; outside 
of it, no explicit absorption is considered. Thus we can 
only follow neutrino emission and the associated neu- 
tronization of matter. However, the up to now micro- 
physically most complete two-dimensional axisymmetric 




Fig. 1. — 3D entropy contours spanning the coordinates planes 
with magnetic field lines (white lines) of the MHD-CCSN simula- 
tion ~ 31 ms after bounce. The 3D domain size is 700 X 700 X 1400 
km. 

study of MHD-CCSN with muhi-group fl ux-limited dif- 
fusion neutrino transport performed by iBurrows et al.l 
(|2007D has shown, that neutrino heating contributes only 
10% 25% to the explosion energy and is therefore sub- 
dominant. This justifies our pragmatic approach at first. 

We employed the pre-coUapse 15M0 model of 
IHeger et al.l ()2005[ ). Although the model provides pro- 
files for rotation and magnetic fields, we use an analytic 
prescription for their distributions and we will comment 
on this choice in Sectional The initial rotation law was 
assumed to be shellular with il{r) — Q,QR^/{r'^ -I- i?g), 
Vto = 11 s~^ and i?o — 1000 km corresponding to an 
initial ratio of rotational energy to gravitational bind- 
ing energy Trot/lW^I = 7.63 x 10~^. For the magnetic 
field we have assumed a homogeneous distribution of a 
purely poloidal field throughout the computational do- 
main of strength 5 x 10^^ G corresponding to an initial 
ratio of magnetic energy to gravitational binding energy 
Tm.J\W\ = 2.63 X 10-8. 

The computed model then undergoes gravitational col- 
lapse and experiences core-bounce due to the stiffening 
of the EoS above nuclear saturation density. Conser- 
vation of angular momentum in combination with the 
collapse leads to a massive spin-up of the core, reach- 
ing Trot / 1 W^ I = 6.81 X 10"^ at bounce, and significant 
rotationally induced deformations. During the collapse 
the magnetic field is amplified by magnetic flux conser- 
vation reaching a central strength of ^ 5 x 10^^ G and 
Tmag/|M^| = 3.02 X 10-4 at bounce. After bounce, differ- 
ential rotation winds up the poloidal field very quickly 
into a very strong toroidal field increasing the mag- 
netic energy/pressure at the expense of rotational energy. 
Consequently, strongly magnetized regions appear near 
the rotational axis with an associated magnetic pres- 
sure quickly reaching and exceeding that of the local gas 
pressure. The Lorentz force then becomes dynamically 
important and matter near the rotational axis is lifted 
from the PNS and drives a bipolar outfiow, i.e., jets are 
launched. The jets rapidly propagate along the rota- 
tional axis and quickly reach the boundary of the initial 
3D domain. In order to follow the jet propagation fur- 
ther, we have continuously extended the 3D domain to a 
final size of 700 x 700 x 1400 km at ~ 31 ms after bounce. 
Figure [T] displays a snapshot at the final time. 

The quickly expanding bipolar jets transport energy 
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Fig. 2. — Time evolution of Ye from a single tracer in green (red) 
for the original data (with neutrino absorption estimate). The blue 
dashed line represents Ye evolved by the network. Vertical dotted 
(dashed) lines correspond to the time when Ue {ve) neutrinosphere 
is crossed; the dot-dashed vertical line when T = 10 GK. 

and neutron rich material outward against the gravi- 
tational attraction of the PNS. We have estimated the 
ejected mass Mej — 6.72 x IQ^^Mq and explosion en- 
ergy i?oxp — 8.45 X 10'*^ erg by summing over the fluid 
cells that are gravitationally unbound. We defined a 
fluid cell as unbound if its total specific energy (inter- 
nal+kinetic-|-magnetic+potential) is positive and if the 
radial velocity is pointing outward. These are admit- 
tedly crude lower bound estimates and these numbers 
were still growing at the end of the simulation. 

3. NUCLEOSYNTHESIS 

The nucleosynthesis calculations are performed with 
a new extended reaction network (jWinteleii 120 lit) 
which represents an advanced (numerically and phys- 
ically) update o f the BasNet network (see, e.g., 
[Thiele mann et al.' (2011)). We us e the reaction 
rates of Rauschcr & Thiclc mannI (|2000l for the FRDM 
mass model). We use the same weak interaction 
rates (electron/positron captures and /?-decays) as in 
lArcones fc Martinez-Pinedd (|2011[ ) . Additionally, we in- 
clude ne utron capture and n eutron induced fission rates 
following IPanov et al.l (120101) and /3-delayed fission prob- 
abilities as described in Pa nov et al.l (j2005l ). 

The tracer particles obtained from the simulation pro- 
vide density, temperature, and electron fraction for the 
nuclear network, as well as position and velocity, from 
the beginning to the end of the simulation [t = tf). Af- 
ter tf , thermodyna mic variable s are ev olved followi ng the 
rescr iption in Fu iimoto et al.l (J2008D . iNishimur a et al.l 
200a) have shown that the details of the expansion 
only have a minor impact on the final abundances. For 
the post-processing we only consider gravitationally un- 
bound tracer particles (see Section [2]). In order to ob- 
tain mass integrated abundances we distribute the total 
ejected mass equally among all ejected tracers. It could 
be shown that this yields very similar results to post- 
processing calculations based o n the conditio ns in the 
unbound cells at the final time ()Winteleiil2011t ). 

The electron fraction is a key input for the nucleosyn- 
thesis and strongly depends on details of the challenging 
neutrino transport. Although neutrino absorption is cru- 
cial to determine the Ye, it is not yet included in the hy- 
drodynamical simulations (where it is expected to have a 
minor impact on the dynamics; see Section[3]). Therefore, 
we present two different nucleosynthesis calculations: (1) 
Fe is taken from the original tracer particles, (2) the ef- 
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Fig. 3. — Ejecta masses vs. Ye for the original simulation without 
neutrino captures (green) and including a simplified prescription 
for neutrino heating (red). The width of a Ye bin is chosen to be 
d,Ye = 0.01. 
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Fig. 4. — Integrated mass fractions for nucleosynthesis calcula- 
tions with (red)and without (green) neutrino heating. Black dots 
represent solar r-process element abundances CSneden et al,.200^ 'l 
scaled to fit the red line at A = 130. 

fects of neutrino absorption on Ye are included also in 
the network. In this second approach, we post-process 
the data from the tracer and use the neutrino informa- 
tion obtained with the leakage scheme (Section [5]). We 
use integrated neutrino luminosities to update the elec- 
tron fraction outside the neutrinosphere (neglecting the 
effect of neutrino energy deposition on the matter tem- 
perature) . The electron fraction of the tracer is evolved 
using approximated rates for the neutrino em ission and 
absorption on nucleons (see, e.g., lJankal ()2001[ )). 

In both of our approaches for Yg , the network calcula- 
tions start when the temperature decreases to 10 GK. 
In the case with post-processing corrections for neu- 
trino absorption, w e also consider neutrino reactions 
from iFrohlich et al.l (J2006D in the nucleosynthesis net- 
work. The required neutrino luminosities and mean en- 
ergies are given by the leakage scheme and assumed to 
be constant and equal to their values at t = tf. 

The evolution of Y^ is presented in Figure [2] for the 
original simulation data and for the estimate of neutrino 
absorption. In the latter, high energy u^ captures on pro- 
tons decrease Y^ outside the i/e-neutrinosphere. However, 
beyond the i^e-neutrinosphere, Ue absorption on neutrons 
dominates and Y^ increases. The fast expansion and the 
relatively low neutrino mean energy limit the effect of the 
absorption. This trend is confirmed also by the network 
for T < 10 GK. 

Figure [3] shows the ejected mass as a function of Yg 
for the original simulation data and for the case in- 
cluding neutrino absorption. These corrections shift the 
peak distribution from ~ 0.15 to ~ 0.17 and broaden 
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it toward higher Y^. In both approaches, at the on- 
set of the nucleosynthesis, the density is stih relatively 
high, p w lO^gcm"^, and the electron fraction rather 

low. Ye ft! 0.15 0.3. This leads to an initial nuclear 

statistical equilibrium composition rich in neutrons and 
neutron-rich nuclei. Such conditions are closer to neu- 
tron star mergers than to the high entropy wind, which 
is characterized by an alpha-rich freeze-out with few nu- 
clei close to stability and small amount of neutrons. 

Mass integrated abundances are presented in Figure |4] 
for the two treatments of the electron frac tion and com- 
■ared to solar r-process abundances (e.g.. iSneden et afl 
20081) ). Nuclei around the second {A = 130) and third 
{A = 195) r-process peaks, as well as the Pb region, 
can be synthesized in the jets. Note that the improve- 
ments in the electron fraction calculation result only in 
minor variations of the abundances. Moreover, the posi- 
tion of the peaks nicely agrees with solar system abun- 
dances. However, the agreement is not equally good 
for all mass numbers as the nuclei between peaks are 
slightly underproduced. Note also that the large trough 

in the mass range 140 160 is mainly due to the strong 

A^ = 82 shell-closure of the FRDM mass model. Simi- 
lar deficiencies were discussed by Nis himura et al.l (|2OO60 
and can be cured by using a different mass model. Nu- 
clei in the region below the second peak are produced 
mainly by charged-particle reactions which occur only 
in the tracer particles with relatively high electron frac- 
tion. This also explains the lower abundances in the 

range A — 80 100 for the case of Y^ without neutrino 

absorption (Figure S]) . This calculation reaches lower Y^ 
which leads to higher neutron densities and a more abun- 
dant third peak and Pb region. 

4. DISCUSSION AND OUTLOOK 

Magneto-rotationally driven supernovae suffer one 
main critique point: the simultaneous presence of fast 
rotation and strong magnetic fields in the progenitor 
bef ore collapse . Rec ent progenitor models computed 
by IHeger et al.l ()2005l ) suggest that these special con- 
ditions are not reached in common evolutionary paths 
of massive stars. However, the low-metallicity models of 
IWoosley & Hcgcr (2006) indicate that for a fraction, per- 
haps 1%, of all massive stars, favorable conditions can ap- 
pear under special circumstances. This rarity of progen- 
itors with these special initial conditions can also be put 
into accordance with the observed scatter of r-process el- 
ement abundances at low metallicity, combined with the 
high r-process production of ^ 6 x 10~^ Mq per event. 
This indicates that these elements h ave been produced i n 
relatively scarce but efficient events (jSneden et al.ll2008l ). 
Therefore, the rare progenitor configurations used here 
could naturally provide a strong r-process site in agree- 
ment with observations of the early galaxy chemical evo- 
lution. 

Currently, the aim to perform long-term global 3D 
simulations of CCSNe, including possibly a sophisti- 
cated radiative transfer of neutrinos, and the aim to 
simulate the local flow instabilities, leading to mag- 



netic field amplification, s eem to be mutually exclusive 
(|Obergaulinger et al.ll2009D . Given this difficulty, we fol- 
lowed here the common approach in the literature by 
taking sufficiently strong initial magnetic fields. By do- 
ing so, it is hypothesized that there is a physical pro- 
cess that can sufficiently quickly amplify the magnetic 
field to dynamic importance by extracting the free rota- 
tional energy in differential rotation. This motivates our 
choice of the high initial poloidal field strength, which by 
flux compression and rotational winding leads to mag- 
netic fields whose magnitude roughly agree with those 
expect ed from the magnetorotatio nal instability at satu- 
ration (|Obergaulinger et al.ll2009D . 

There are additional shortcomings in the treatment 
presented here. This includes that thermodynamic prop- 
erties of tracer particles are only extrapolated beyond 
the end of the MHD calculation and the nucleosyn- 
thesis results were not yet tested with several nuclear 
mass models. But the main outcome of this investiga- 
tion is that full 3D calculations can support the emer- 
gence of bipolar jets and that these are not artifacts 
of up to now axisymmetric approaches (jNishimura et all 
,2^006; Fujimoto ct al. 2008). Such explosions, resulting 
from the individual evolution of massive stars rather 
than complicated binary histories of neutron star merg- 
ers, could explain strong r-process features during early 
galactic evolution as observed in many low metallicity 
stars. The r-process production is very efficient, a fac- 
tor of 100 more than expected on average from super- 
novae, if they would be responsible for solar r-process 
abundances. The question remains: whether magnetars 
— with the magnetic fields required for this outcome — 
result from about 1% of all core collapses or are rarer 
event s. Ap parently present knowledge permits this op- 
tion (jKouv eliotou ct al. 1998; Kramer 2009). 

The magnetorotationally driven supernova simulation 
presented here provides a scenario for a strong r-process 
site seemingly consistent with observations of the early 
chemical evolution of our Galaxy. Future studies should 
investigate the effects of the progenitor configuration 
(like mass, magnetic field and rotation rate), the nuclear 
input (in particular the mass formula), and a more con- 
sistent treatment of neutrino interactions on the robust- 
ness of the conclusions drawn here. 

We thank T. Rauscher and I. Panov for useful discus- 
sions and help with reaction rates and fission. We ac- 
knowledge support and funding by the Swiss National 
Science Foundation (SNSF). A. A. was supported by 
Feodor Lynen Fellowship (Humboldt Foundation) and 
by the Helmholtz-University Young Investigator grant 
no. VH-NG-825. The authors are additionally supported 
by EuroGENESIS, a collaborative research program of 
the European Science Foundation (ESF), CompStar, a 
research networking program of the ESF, and the FP7 
ENSAR project, funded by the European Commission. 
This work was supported by a grant from the Swiss Na- 
tional Supercomputing Centre (CSCS) under project ID 
s332. 



REFERENCES 



Arcones, A., Janka, H.-T., & Scheck, L. 2007, A&A, 467, 1227 



Arcones, A. & Martinez-Pinedo, G. 2011, Phys. Rev. C, 83, 
045809 



r-process in 3D MHD Jets 



Argast, D., Samland, M., Thielemann, F.-K., & Qian, Y.-Z. 2004, 

A&A, 416, 997 
Bethe, H. A. 1990, Rev. Mod. Phys., 62, 801 
Bisnovatyi-Kogan, G. S., Popov, Y. P., & Saniochin, A. A. 1976, 

Ap&SS, 41, 321 
Bruenn, S. W. 1985, ApJS, 58, 771 
Burrows, A., Dessart, L., Livne, E., Ott, C. D., & Murphy, J. 

2007, ApJ, 664, 416 
Burrows, A., Nordhaus, J., Almgren, A., & Bell, J. 2011, 

Comput. Phys. Commun., 182, 1764 
Cameron, A. G. W. 2003, ApJ, 587, 327 

Farouqi, K., Kratz, K.-L., Pfeiffer, B., et al. 2010, ApJ, 712, 1359 
Fischer, T., Whitehouse, S. C., Mezzacappa, A., Thielemann, 

F.-K., & Liebendorfer, M. 2010, A&A, 517, ABO 
Frohlich, C., Martinez-Pinedo, G., Liebendorfer, M., et al. 2006, 

Phys. Rev. Lett., 96, 142502 
Fujimoto, S., Nishimura, N., & Hashimoto, M. 2008, ApJ, 680, 

1350 
Heger, A., Woosley, S. E., & Spruit, H. C. 2005, ApJ, 626, 350 
Hoffman, R. D., Woosley, S. E., & Qian, Y.-Z. 1997, ApJ, 482, 

951 
Hiidepohl, L., Miiller, B., Janka, H.-T., Marek, A., & Raffelt, 

G. G. 2010, Phys. Rev. Lett., 104, 251101 
Janka, H.-T. 2001, A&A, 368, 527 
Kappeli, R., Whitehouse, S. C., Scheidegger, S., Pen, U.-L., & 

Liebendorfer, M. 2011, ApJS, 195, 20 
Kouveliotou, C., Dieters, S., Strohmayer, T., et al. 1998, Nature. 

393, 235 
Kramer, M. 2009, in lAU Symp. 259, ed. K. G. Strassmeier, A. G. 

Kosovichev, & J. E. Beckman (Cambridge: Cambridge Univ. 

Press), 485 
Lattimer, J. M. & Swesty, F. D. 1991, Nuc. Phys. A, 535, 331 
LeBlanc, J. M. & Wilson, J. R. 1970, ApJ, 161, 541 
Lentz, E. J., Mezzacappa, A., Bronson Messer, O. E., et al. 2012, 

ApJ, 747, 73 
Liebendorfer, M., Pen, U., & Thompson, C. 2005, Nucl. Phys. A, 

758, 59 
Liebendorfer, M., Rosswog, S., & Thielemann, F.-K. 2002, ApJS, 

141, 229 
Liebendorfer, M., Fischer, T., Hempel, M., et al. 2010, Prog. 

Theor. Phys. Suppl., 186, 87 



Marek, A., Dimmelmeier, H., Janka, H.-T., Miiller, E., & Buras, 

R. 2006, A&A, 445, 273 
Meier, D. L., Epstein, R. I., Arnett, W. D., & Schramm, D. N. 

1976, ApJ, 204, 869 
Meyer, B. S. & Brown, J. S. 1997, ApJS, 112, 199 
Moiseenko, S. G. & Bisnovatyi-Kogan, G. S. 2008, Int. J. Mod. 

Phys. D, 17, 1411 
Miiller, B., Janka, H.-T., & Dimmelmeier, H. 2010, ApJS, 189, 

104 
Miiller, E. & Hillebrandt, W. 1979, A&A, 80, 147 
Nishimura, N., Fischer, T., Thielemann, F.-K., et al. 2011, 

arXiv:1112.5684 
Nishimura, S., Kotake, K., Hashimoto, M., et al. 2006, ApJ, 642, 

410 
Obergaulinger, M., Cerda-Duran, P., Miiller, E., & Aloy, M. A. 

2009, A&A, 498, 241 
Panov, I. v., Kolbe, E., Pfeiffer, B., et al. 2005, Nucl. Phys. A, 

747, 633 
Panov, I. v., Korneev, I. Y., Rauscher, T., et al. 2010, A&A, 513, 

A61 
Pruet, J., Hoffman, R. D., Woosley, S. E., Janka, H.-T., & Buras, 

R. 2006, ApJ, 644, 1028 
Rauscher, T. & Thielemann, F.-K. 2000, At. Data Nucl. Data 

Tables, 75, 1 
Rosswog, S. & Liebendorfer, M. 2003, MNRAS, 342, 673 
Sagert, L, Fischer, T., Hempel, M., et al. 2009, Phys. Rev. Lett., 

102, 081101 
Scheidegger, S., Kappeli, R., Whitehouse, S. C, Fischer, T., & 

Liebendorfer, M. 2010, A&A, 514, A51-I- 
Sneden, C, Cowan, J. J., & Gallino, R. 2008, ARA&A, 46, 241 
Takiwaki, T., Kotake, K., & Sato, K. 2009, ApJ, 691, 1360 
Thielemann, F.-K., Arcones, A., Kappeli, R., et al. 2011, Prog. 

Part. Nucl. Phys., 66, 346 
Thompson, T. A., Burrows, A., & Meyer, B. S. 2001, ApJ, 562, 

887 
Wanajo, S. 2006, ApJ, 650, L79 

Wanajo, S., Janka, H.-T., & Miiller, B. 2011, ApJ, 726, L15 
Winteler, C. 2011, PhD thesis. University of Basel 
Woosley, S. E. & Heger, A. 2006, ApJ, 637, 914 



